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ABSTRACT 

A few rare halo giants in the range [Fe/H] ~ —2.9 ±0.3 exhibit r-process element abundances 
that vary as a group by factors up to [r/Fe] ~ 80, relative to those of the iron peak and below. 
Yet, the astrophysical production site of these r-process elements remains unclear. We report 
initial results from four years of monitoring the radial velocities of 17 r-process-enhanced metal- 
poor giants to detect and characterise binaries in this sample. We find three (possibly four) 
spectroscopic binaries with orbital periods and eccentricities that are indistinguishable from those 
of Population I binaries with giant primaries, and which exhibit no signs that the secondary 
components have passed through the AGB stage of evolution or exploded as supernovae. The 
other 14 stars in our sample appear to be single - including the prototypical r-process-element 
enhanced star CS 22892-052, which is also enhanced in carbon, but not in s-process elements. 
We conclude that the r-process (and potentially carbon) enhancement of these stars was not a 
local event due to mass transfer or winds from a binary companion, but was imprinted on the 
natal molecular clouds of these (single and binary) stars by an external source. These stars are 
thus spectacular chemical tracers of the inhomogeneous nature of the early Galactic halo system. 

Subject headings: Galaxy: halo — Stars: Population II — binaries: spectroscopic — techniques: radial 
velocities 

1. INTRODUCTION 
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history of the early Galaxy. Carbon (and often 
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for stars with [Fe/HJ < - 2.5 (|Beers k Christlieb 



2005; ICarollo et al.N2011l ). but the overall abun- 



dance ratios of elements up to the iron peak are 
well esta blished by [Fe/H] ~ —2.5, with very sma ll 



scatter (|Cavrel et all 12004 lArnone et"al1 [§005) 



The full range of r-process elements was also 
in place at this stage in the chem ical evolution 
of th e Galaxy (see the review by ISneden et al 



2008). Other neutron-capture processes, such 



as the s-process, began to contribute significant 
amo unts of heavy element s from [Fe/H] < —2.5 
(e.g, Simmerer et al . 20040. Hence, the processes 
by which the chemical elements were produced and 
recycled into the Galactic halo system are, at least 
in a broad-brush sense, well understood. 

However, in a small fraction of stars in the 
abundance range —3.2 < [Fe/H] < —2.6, the uni- 
form abundance pattern of the r-process elements 
as a group is enhanced by factors up to ~ 80 
relative to that of the iron-peak and lighter ele- 
ments. Spectroscopic analyses of such stars with 
8m-class telescopes have provided precise and de- 
tailed abundances of many r-proces s elements as 



a key to understanding t heir origin (|Sneden et al 



2008; ICowan et alj|201ll) . Yet, the likely produc 



tion site(s) of the r-process elements, as well as 
the mechanisms by which their abundances rela- 
tive to the "standard" halo composition could vary 
so strongly from star to star in the early Galaxy, 
remain essentially unknown. 

Explanations of this diversity fall into two main 
classes: Inhomogeneous enrichment and incom- 
plete mixing of the ISM by the first generation(s) 
of stars, or later, local pollution of neutron-capture 
elements by a b inary companion of t he pr esently 
observed star ( Qian fc Wasserburel 200 ll) . In 
the first case, the very uniform abundance pat- 
tern of the a- and iron-peak elements is diffi- 
cult to reconcile with the predictions of models 
of stochastic star formation and enrichment in 



the early Galaxy (jArnone et al.l 120051 ) . This con- 
flict would be resolved in the second case, but 
unlike the s-process-element enhanced Ba and CH 
gian ts, which are known to all be long-period bina- 
ries jMcClure fc WoodsworthT lQiffl: iJorissen et al 



19981 ). this conjecture is so far without observa- 



^ote that Rocdcrcr ct al. have argued that the ma- 

jority of s— process element enhancement may have been 
delayed until [ Fe/H] ~ —1.4. See also the discussion in 
iBisterzo et"aD (20 111 ). 



tional foundation. 

Here we report the first results from four years 
of precise radial-velocity monitoring, performed in 
order to assess the binary frequency of a sample 
of 17 r-process-element enhanced VMP and EMP 
giants. Our results provide strong new constraints 
on the nature of the r-process production site(s) 
and on the use of these stars as tracers of the 
star formation and/or merger history of the early 
Galaxy. 

2. Sample Definition and Observations 

Our program stars were drawn from the HERES 

search fo r r-process-element e nhanced stars (jChristlieb et al 



2004; iBarklem et al.l 120051 ) . supplemented by 



earlier and l ater d iscoveries as summarised by 
Havek et al.l (|2009l ). Only stars north of decli- 
nation ~ —25° and brighter than V ~ 16.0 are 
accessible for study with the Nordic Optical Tele- 
scope (NOT) on La Palma, resulting in a total 
sample of 17 stars. Eight of these are in the r- 
I class (+0.3 < [Eu/Fe] < +1.0), as defined by 
Beers fc Christliebl (|2005l ). and 9 are in the r-II 
class ([Eu/Fe] > +1.0). Table Q] lists the pro- 
gram stars, in order of increasing [r/Fe] ratio, in 
order to highlight the continuum of r-process en- 
hancements that exists; the division into r-I and 
r-II classes appears to be merely one of conve- 
nience. Stars with measured U abundances, and 
those found here to be spectroscopic binaries, are 
indicated in the Table. 

As the chemical compositions of the targets 
are already known - some in great detail - our 
observations were designed to just yield precise 
radial velocities as efficiently as possible, over 
a time span of several years. To this end, we 
obtained high-resolution spectra (R ~ 45, 000) 
with a S/N ratio ~ 10, using the bench- mounted, 
fibre-fed echelle spectrogra ph FIES at the NOT 
(|Djupvik fc Andersenll2010l ). which is installed in a 
separate, temperature-controlled underground en- 
closure. The useful wavelength range covered by 
these spectra is 4000 - 7000 A. 

Our goal was to reach an accuracy of 100 — 



2 Note that the r-II and r-I stars do occupy different, but 
overlapping, regions of metallicity space. The r-II stars are 
found only at very low metallicity, while the r-I stars extend 
over the ra nge —3.0 < [Fe/H] < —0.5; see the bottom panel 
of Fig. 3 in lAoki et atl (2010T ). 
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200 m s" 1 per observation, except perhaps for the 
faintest program stars. The radial-velocity zero 
point was checked with standard stars on every 
night of observation, and found to be reproducible 
to better than 45 m s _1 per observation over a 
five-year period. Thus, the accuracy of the radial 
velocities of the program stars is not limited by 
the instrument. 

Our initial assu mption, by analogy with the 
Ba and CH giants (|Jorissen et al.lll998l ). was that 
any spectroscopic binaries in the sample would 
likely have orbits of long period, low eccentric- 
ity, and small amplitude - i.e., small, slow veloc- 
ity variations. Thus, our strategy was to observe 
these stars at roughly monthly intervals, weather 
permitting, and then adapt the frequency of the 
observations to follow any objects with radial- 
velocity variations detected in the initial data. 
Observations began in April 2007, and were con- 
tinued on 51 nights through September 2011, for 
a total of ~ 234 spectra, an average of 14 per 
star. Faint stars at far southern declinations re- 
quire ideal conditions, and were observed less fre- 
quently than average. 

3. Data Reduction and Analysis 

The entire set of reductions of the raw spec- 
tra (bias subtraction, division by a flat-field ex- 
posure, cosmic ray removal, 2-D order extraction, 
and wavelength calibration) was performed with a 
program develo ped and ex t ensive ly tested on exo- 
planet hosts by iBuchhavd (l2010h . For the fainter 
stars, it was found preferable to divide the long ex- 
posures into three pieces, and remove cosmic ray 
hits by median filtering. 

Radial velocities were then derived from the re- 
duced spectra by a multi-order cross-correlation 
procedure. This operation is the most difficult 
step in the analysis because (i) these stars are ex- 
tremely metal-poor and chemically peculiar; and 
(ii) the individual spectra have low S/N ratios. 
Thus, selecting an optimum template spectrum for 
each star is no trivial task. Noting that the pri- 
mary objective of the analysis is to measure small 
velocity variations rather than absolute values, we 
have experimented with three types of template 
spectra: (a) the highest S/N spectrum of each 
star; (b) the velocity-shifted and co-added mean 
spectrum of each star, and (c) a synthetic spec- 



trum consisting of 5 functions at the (solar) wave- 
lengths of the strongest visible lines. The choice 
of template for each star was then guided by the 
consistency of the resulting velocities. Templates 
(b) and (c) were generally found to give the most 
consistent results, the latter also yielding veloci- 
ties on a reliable absolute scale. 

In summary, our final procedure yielded radial 
velocities with a standard deviation of < 100 m 
s _1 for the brighter and more metal-rich stars, ris- 
ing to ~ 300 — 1000 m s^ 1 for the faintest stars 
with the weakest spectral features. 

4. Binary Detection and Orbit Determi- 
nation 

Fourteen of our stars exhibit no significant vari- 
ation in radial velocity over the period covered by 
our observations, inclu ding any earlier vel ocities 
reported from HERES (jBarklem et al.ll2005l ). Lin- 
ear and parabolic fits of the run of velocities vs. 
time were made to check for any long-term trends, 
but they were generally negligible within the un- 
certainties; a few of the stars will be kept under 
continued surveillance. 



Th e star HE 0442-1234 was shown bv lBonifacio et al 
(2010) to be a probable long-period spectroscopic 
binary. Our new data enabled us to complete 
the orbit of this star, as well as for the newly- 
discovered binaries HE 1044-2509 and HE 1523- 
0901. The orbital elements for these stars are 
listed in Table [3J and the observed and computed 
velocity curves are shown in Fig. [T] 

In general, the individual component masses 
cannot be derived directly, but assuming a stan- 
dard value of 0.8 M Q for the mass of a halo giant 
allows us to estimate a minimum mass for the un- 
seen companion. For HE 0442-1234, this is 0.67 
M & if the inclination i ~ 90 deg, and it cannot be 
much larger for the secondary star to remain in- 
visible in the spectrum. In the other two systems, 
a secondary star on the main sequence could be 
of similar or lower mass. Assuming M% = 0.6M Q 
leads to the estimates of i given in Table [2] note 
that the orbit of HE 1523-0901 is seen nearly face- 
on. These, in turn, lead to estimates of the size 
(volume equivalent radius) of the Roche lobes of 
the unseen stars, which are not very sensitive to 
the adopted geometry. 
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Table 1 

Stars Monitored for Radial Velocity Variation 



Star 


RA (J2000) 


Doc (.T2000) 


V 


B - V 


[Fe/H] 


[r/Fe] 


Nobs 


Remarks 


HE 0524-2055 


05:27:04 


-20:52:42 


14.01 


0.87 


-2 


.58 


+0.49 


8 


r-I 


HE 0442-1234 


04:44:52 


-12:28:46 


12.91 


1.07 


-2 


.41 


+0.52 


23 


r-I, SB 


HE 1430+0053 


14:33:17 


+00:40:49 


13.69 


0.58 


-3 


.03 


+0.72 


19 


r-I 


CS 30315-029 


23:34:27 


-26:42:19 


13.66 


0.91 


-3 


.33 


+0.72 


9 


r-I 


HD 20 


00:05:15 


-27:16:18 


9.07 


0.54 


-1 


.58 


+0.80 


9 


r-I 


HD 221170 


23:29:29 


+30:25:57 


7.71 


1.02 


-2 


.14 


+0.85 


23 


r-I 


HE 1044-2509 


10:47:16 


-25:25:17 


14.34 


0.66 


-2 


.89 


+0.94 


13 


r-I, SB 


HE 2244-1503 


22:47:26 


-14:47:30 


15.30 


0.60 


-2 


.88 


+0.95 


12 


r-I 


HE 2224+0143 


22:27:23 


+01:58:33 


13.68 


0.71 


-2 


.58 


+ 1.05 


18 


r-II 


HE 1127-1143 


11:29:51 


-12:00:13 


15.88 




-2 


.73 


+ 1.08 


12 


r-II 


HE 0432-0923 


04:34:26 


-09:16:50 


15.16 


0.73 


-3 


.19 


+ 1.25 


14 


r-II 


HE 1219-0312 


12:21:34 


-03:28:40 


15.94 


0.64 


-2 


.81 


+ 1.41 


6 


r-II 


CS 22892-052 


22:17:01 


-16:39:26 


13.21 


0.80 


-2 


.95 


+ 1.54 


17 


r-II 


CS 29497-004 


00:28:07 


-26:03:03 


14.03 


0.70 


-2 


.81 


+ 1.62 


8 


r-II 


CS 31082-001 


01:29:31 


-16:00:48 


11.66 


0.76 


-2 


.78 


+ 1.66 


17 


r-II, U 


HE 1523-0901 


15:26:01 


-09:11:38 


11.10 


1.10 


-2 


.95 


+ 1.80 


17 


r-II, U, SB 


HE 1105+0027 


11:07:49 


+00:11:38 


15.64 


0.39 


-2 


.42 


+ 1.81 


9 


r-II 



Note. — Remarks: U indicates that uranium has been detected; SB indicates a confirmed spectroscopic binary. 



Table 2 

Orbital Elements for the Detected Binary Stars 



Element HE 0442-1234 HE 1523-0901 HE 1044-2509 



P (d) 2513.38+4.46 302.78+0.78 36.57+0.20 

K(kms- 1 ) 12.50+0.20 0.399+0.008 27.48+0.22 

e 0.760+0.001 0.23+0.129 0.000+0.000 

7 (kms" 1 ) 236.16+0.20 162.50+0.10 359.28+0.55 

f(m) (M ) 0.1396+0.0012 1.59 10~ 6 ±0.54 10~ 6 0.0823+0.0027 

asini (R e ) 403.2+1.2 2.21+0.17 20.24+0.24 

i (for M 2 ~ 0.6M Q ) 88 1.5: 65 

RRoche (Re) 65° 48 13 

a (1 obs., km s" 1 ) 0.28 0.11 0.98 



a Size at periastron. 
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5. Discussion 




Fig. 1. — Observed and computed spectroscopic 
orbits for HE 0442-1234 (top), HE 1044-2509 
(middle) and HE 1523-0901 (bottom). Orbital el- 
ements are listed in Table [5] Blue dots: FIES 
velocities; red dots: Earlier data. 



Our results enable us to address several issues of 
importance for understanding the likely astrophys- 
ical site(s) of the r-process, and also shed some 
light on the early Galactic production of carbon, 
as described below. 

5.1. Binary Frequency 

With three binaries detected in a sample of 
17 stars, the binary frequency of the r-process- 
enhanced stars is 18%. The star HE 2327-5642, 
an r-II star showing a pos sibly variable veloc- 



ity ( Mashonkina et al. 2010h . is below the south- 



ern limit of our sample, but is another candi- 
date. This is fully consistent with the ~20% bi- 
nary fr equency determine d for normal cluster gi- 
ants bv IMermilliodl (jl996l) . An additional 1-2 fu- 
ture discoveries in a larger sample might boost 
the frequency to perhaps as much as 25%, but 
there is clearly no evidence that all r-process- 
ele ment enhanced stars are bin aries, as speculated 
bv lQian fc Wasserbure ( 2001 ) and others. 

Within the limitations imposed by the small 
sample, there also seems to be no difference in 
the binary population of the r-I and r-II classes. 
Similarly, of the two stars with measured U abun- 
dances, CS 31082-001 exhibits the so-called "ac- 
tinide boost" (an over-abundance of Th and U 
relative to third-peak r-process elements such as 
Eu) and is a single star, while HE 1523-0901 is a 
binary and shows no actinide boost. Remarkably, 
the C-enhanced prototypical r-II star CS 22892- 
052 also seems no t to be a binary, de s pite th e ear- 
lier suggestion bv IPreston fc Snedenl ( 2001 ). indi- 
cating that its C content was not produced in a 
former AGB companion. Thus, membership in a 
binary system appears to be decoupled from de- 
tails associated with these particular abundance 
variations. 

5.2. Orbital Properties 

The periods and eccentricities of our three con- 
firmed binaries are entirely consistent with those 
of the sample of chemically norma l Population I 
giant binaries by IMermilliodl (|1996) : The longest- 
period orbit is highly eccentric, and the shortest- 
period orbit is well below the limit of ~ lOOd for 
tidal circularization. Note that the old, metal- 
poor CH stars typically have circular orbits for 
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periods up to ~ lOOOd. Moreover, the secondary 
Roche lobes are too small to accommodate typ- 
ical AGB stars of ~ 200 Rq, if the secondary 
stars passed through this phase of their evolution. 
This is again consistent with the lack of s-process- 
element enhancements in these stars. 

5.3. The Astrophysical Site of the r- 
Process 

Proving that the r-process-element production 
did not occur in binaries provides no direct evi- 
dence of the nuclear physics process at the produc- 
tion site above that afforded by the detailed abun- 
dance analyses of these stars. Models attempt- 
ing to explain the observed r-process abundance 
patterns fall into two classes, core-collapse super- 
novae (S Ne) and merging binary neutron star s 
(see, e.g.. lArgast et all 12004 iGoriely et aT]|201ll) 
also d iscussed comprehensively by ISneden et al.l 
(120081). and in the context of CS 31082-001 by 
Barbuv et abl (j201ll ). 

The key feature to be highlighted here is that 
the newly-formed r-process elements were added 
in variable, but internally consistent proportions, 
in the otherwise constant chemical composition 
of the next generations of EMP and VMP stars. 
That the r-process-element abundances varied so 
strongly in the clouds from which these stars 
formed indicates that the r-process elements were 
not simply uniformly dispersed, together with all 
lighter elements, in the supernova explosions that 
are a common feature of all astrophysical models 
for the r-process. Ejection in a jet or beam directly 
from the nascent neutron star(s) seems the most 
natural scenario for achieving this. The chemical 
composition of the progenitor, and the varying dis- 
tance and direction of the jet from the next cloud, 
would then explain, in a natural way, the continu- 
ously varying proportions of r-process elements to 
the bulk composition of the following generation 
of EMP stars. 

5.4. r-Process-Enhanced Stars as Chemi- 
cal Tags of the Early Galactic ISM 

In this scenario, most stars would receive a 
"standard" dose of r-process elements; stars with 
r-process-element abundances exceeding the aver- 
age by a factor of > 2—3 would be seen as r-process 
rich, while those below the average (exemplified 



by H P 122653, IWestin et all l2000t iHonda et al 
2006) would have been bypassed by the r-process 
ejection and appear as r-process poor. This lat- 
ter group may be as numerous as the former, 
but without spectacularly strong spectral lines to 
call attention to them. The recent discovery by 



Aoki et all (|2010l ) of a cool, EMP mam-sequence 
dwarf with highly r-process-enhanced elemental 
abundance ratios, consistent with classification as 
an r-II star, would obviate any model to explain 
r-process enhancements as only due to some at- 
mospheric anomaly confined to red giants. 

The existence of r-process-element enhanced 
(or depleted) stars in a narrow range of metal- 
licity near [Fe/H] ~ —3 would imply that such 
anisotropic SN explosions only appeared at a 
certain "chemical time", and that the ISM was 
quickly fully mixed soon thereafter. The spectac- 
ular abundance anomalies of these stars can thus 
be used as extreme examples of "chemical tags" 
of the sites and times of their formation. The 
r-I/r-II classification is just a coarse tool to indi- 
cate the degree of enhancement. However, as r-I 
stars are also found at higher metallicity than the 
r-II stars, they may have formed from clouds that 
were further enriched by material of "normal" halo 
composition. 

The different processes responsible for the 
light and heavy r-process elements ( see, e.g., 



Francois et al.l 120071 : iMontes et al.l 120071 ). as well 



as the existence of stars with and without an ac- 
tinide boost, remain to be explained by further 
modeling, but the binary properties of EMP stars 
apparently played no role in this context: Binaries 
formed with similar p roperties as in chemically 
norm al stars 
20081) . 



(e.g., iGonzalez Hernandez et al 



Finally, it is remarkable that the prototypical 
r-II star CS 22892-052 is single and significantly 
enhanced in carbon, which has been assumed to 
originate in long-period binary stars together with 
the s-process elements, which are not observed in 
CS 22892-052. This casts doubt on the accepted 
explanation for the synthesis of C in the early 
Galaxy as due primarily to pollution by former 
AGB binary companions, and suggests the synthe- 
sis of C, N, and O in earlier, rapidly rotating mas- 
sive stars as one a t tracti ve alternative (see, e.g., 
Mevnet et aDl2006ll201oh . 
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6. Conclusions 

Eighty percent of our program r-I and r-II stars 
exhibit no detectable radial-velocity variations, 
while three stars are binaries with well-determined 
orbits (Table [2]), typical of systems with giant pri- 
maries, but no AGB secondary stars. Thus, the 
binary population among these stars is normal, 
and binary stars play no special role in producing 
the r-process elements and injecting them into the 
early ISM. The case of CS 22892-052 suggests that 
this may be true for the early synthesis of carbon 
as well. 

We conclude that whatever progenitors pro- 
duced the r-process elements (and carbon) were 
extrinsic to the EMP and VMP stars we observe 
today. These elements were likely ejected in a col- 
limated manner, and make these stars archetypical 
chemical indicators of their formation sites in the 
early Galaxy. 

This paper is based on observations made with 
the Nordic Optical Telescope, operated on the is- 
land of La Palma jointly by Denmark, Finland, 
Iceland, Norway, and Sweden, in the Spanish Ob- 
servatorio del Roque de los Muchachos of the In- 
stitute de Astrofisica de Canarias. We thank 
Drs. Piercarlo Bonifacio, Luca Sbordone, Lorenzo 
Monaco, and Jonay Gonzalez Hernandez for alert- 
ing us to the binary nature of HE 0442-1234 and 
for allowing us to include their velocities in our 
orbital solution, Dr. G. Torres for computing the 
orbit of HE1523-091, and Dr. Paul Barklem for 
providing the observing dates of the HERES spec- 
tra. We also thank several NOT students for ob- 
taining most of the NOT observations for us in 
service mode. J. A. and B.N. acknowledge sup- 
port from the Danish Natural Science Research 
Council, and L.A.B from the Carlsberg Founda- 
tion. T.C.B. acknowledges partial funding of this 
work from grants PHY 02-16783 and PHY 08- 
22648: Physics Frontier Center/ Joint Institute for 
Nuclear Astrophysics (JIN A), awarded by the U.S. 
National Science Foundation. 
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